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Abstract. By investigating the correlations between dust column density as inferred 
from infrared data and the observed colours of celestial objects at cosmological 
distances with small colour dispersion, we constrain the properties of Milky Way 
dust. Results derived using colours of quasars, brightest central galaxies and large 
red galaxies are broadly consistent, indicating a proportionality constant between 
the reddening E{B — V) = Ab — Ay and the dust column density (given 
in units of MJy/sr) oi p — E{B — V)/D'^ ~ 0.02 and a reddening parameter 
Rv = Av I E{B — V) ^ i with fractional uncertainties of the order 10%. The data 
does not provide any evidence for spatial variations in the dust properties, except 
for a possible hint of scatter in the dust extinction properties at the longest optical 
wavelengths. 



1. Introduction 

Astronomical data need to be corrected for extinction from dust in the Milky Way. 
Typical (albeit with large variations) values for the reddening E{B — V) = Ab — Ay 
is 0.02 with Ay ~ 0.06, see Fig. [H [|| In this paper, we reinvestigate the relation 
between the Milky Way reddening and the dust maps generated in Schlegel et al. (l2| . 
using colours of quasars (QSOs) observed within the Sloan Digital Sky Survey (SDSS). 
We also compare with results from SDSS observations of Brightest Central Galaxies 
(BCGs) and Luminous Red Galaxies (LRGs). Compared with the original investigation 



in Schlegel et al. [12| , the improved statistics and multiple band photometry potentially 
allow us to decrease the errors on a linear relation while simultaneously constraining 
the properties of Milky Way dust, in this case the total to selective extinction ratio 
Rv = Av/E{B~V). 

We also perform a simple check to investigate if there are deviations from the linear 
relation in that regions with high dust column density have more or less reddening than 
the linear prediction. No such effect is found. Neither do we find any firm evidence of 
scatter in the dust properties as a function of the line-of-sight. 

I The dust maps of Schlegel et al. [lij with degree resolution peak at E{B — U) ~ 0.016 mag, with a 
mean value of 0.048 mag, and a median value of 0.028 mag. 
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Figure 1. Galactic dust extinction represented as — V^) with degree resolution for 
galactic latitudes |&| > 20°, together with the corresponding probability distribution 
function. Data taken from Schlegel et al. [l^ . 



While noting that there is still room for non-linear behaviour as well as scatter in 
dust properties at angular positions and scales not probed by the data, we conclude that 
at the positions and scales studied in this paper, it is a fair assumption that the reddening 
is linear to the dust column density (in units of MJy/sr) with proportionality 
constant p = E{B - V)/D^ ~ 0.02 and Ry ~ 3. The relative errors are of the order 
of 10 % and are mainly limited by uncertainties in the filter response functions and the 
source spectral templates. 

This agrees well with the results of Berry et al. [l| where SDSS photometry of 73 
million stars is used to obtain the dust extinction towards each star giving a three- 
dimensional dust map of the Milky Way in the direction of the observed stripes. The 
results from this analysis are consistent with Ry = 3.0, with a statistical and systematic 
error of 0.1 and a scatter of 0.2. 



In Peek and Graves 10|, the colours of passively evolving galaxies from the SDSS 



are used to obtain corrections to the results of Schlegel et al. [12[. They find that 
the reddening in regions of low dust temperature are underpredicted with maximum 
deviations in E{B — V) of 45 millimagnitudes. 

In the following, all limits are given as 95.4% confidence levels (CL) corresponding 
to = Xmin + 3-98 for one parameter and = xLin + 6-16 for two parameters. 

2. Method 



We use the temperature corrected dust maps from Schlegel et al. [12| to derive a 
reddening map by correlating the dust maps to the observed colours of celestial objects. 
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preferably with a small intrinsic colour dispersion. The dust map, given as the dust 
column density in units of MJy/sr, is derived from COBE/DIRBE data allowing 
for a resolution of approximately one degree. We note that in Kohyama et al. [6i], IRAS 
data is used to reach a resolution of five arc minutes in the Cygnus region. However, no 
all sky map data is available for these data. 

As a first attempt, we assume that the observed colour, e.g. E{B — V), is 
proportional to the dust column density 

E{B - V) = pD'^ + m, (1) 

where we set out to derive the proportionality constant p, and m represents an arbitrary 
offset that is marginalised over (possibly due to, e.g., calibration errors giving the 
wrong mean colour) and the dust column density is given in units of MJy/sr, making 



p dimensionless. In Schlegel et al. [12[, a value of the reddening parameter Ry = 3.1 
was assumed, yielding a value of the proportionality constant p = 0.0184 ± 0.0014 for 
Mg2-cahbrated BCGs. 

Note that this can in principle be done for any set of filters, X and Y, 

E{X-Y) = kxYD'^ + mxY, (2) 
and that multiple colour data thus can be used to infer dust properties. 

2.1. Fitting for Ry and p 

Since SDSS data allow for multiple colour combinations, we are able to fit also 
for the properties of the Milky Way dust and relax any prior assumptions on the 
reddening parameter Ry. The SDSS [u,g,r,i,z] filters have central wavelengths at 
[3543,4770,6231,7625,9134] A and we thus expect g — r and/or m — r to correspond 
closest to B — V. 

For each source and filter, we calculate the transmitted flux using a redshifted 
source template. This is done for a grid of dust models in order to compute how each 
dust models will affect a given colour. Note that in doing this, we need to calculate 
the transmitted flux for each dust model separately since /^-corrections will depend on 
the dust model. In practice, we have created tables that list the necessary fc-corrections 
over a grid of redshifts. Ay and Ry for each filter that we interpolate between. 

After calculating how a given dust model affects each source template in each filter, 
we can compare with the actual observed colours, E{X — Y), allowing us to constrain 
the parameters p and Ry. 

Normalizing such that zero column density corresponds to zero reddening, i.e. 
m = in Eq. [H defining /(A) = A{\)/Ay and using 

Ay = pRyD^ , (3) 



we can write 

Al- A2 _ Am 



pRvifi - f2) , (4) 
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and 

(Am \ 
—^j = dpRvifi - /2) + pdRvifi - f2) + pRvdUi - f2) , (5) 

showing how uncertainties in , Am and the dust extinction curve /(A) translates into 
uncertainties in p and Ry- In the following, we will make use of the Milky Way extinction 
curves /(A) as parametrised in Cardelli et al. 0] (CCM) and Fitzpatrick ^3] (FTZ), 
plotted in Fig. [3] together with the (arbitrary normalised) SDSS filter transmission 
curves. 

We now turn to the expected degeneracy between Ry and p and how this can be 
broken using multiple colour data. In the CCM-model, we have / = a + b/Ry and 
fi- f2 = (ai - 02) + (fei - b2)/Rv = Aa + Ab/Ry. Setting d (Am/D^) = in Eq. E 
and putting d{Aa) = d{Ab) = 0, i.e. no systematic error in the extinction curves, we 
obtain 

dp _ -dRy 

p ~ Rv + Ab/Aa ■ ^ ^ 

This can be integrated to give 

i?y,o + Ab/Aa 
Ry + Ab/Aa 

describing the relation between different values of Ry and p giving the same extinction 
per unit column density Am/D'^. The degeneracy lines for the SDSS standard colours 
u — g (dark red), g — r (red), r — i (orange) and i — z (yellow) are plotted in Fig. Ej 
displaying the power of multiple colour data to break the degeneracy between p and Ry. 
We see that p is anti-correlated with Ry, if the difference {Ax — Ay) /Ay is decreasing 
slower than 1/Ry with Ry, which is the case for the three reddest SDSS standard 
colours. The slope of these curves compare well with the confidence contours derived 
from data in the forthcoming sections. 



3. Error analysis 

We first note that constant filter calibration errors do no affect our results, since we are 
marginalising over any constant offset between the observed and theoretically calculated 
colours. However, we may get systematic effects if the spectral template is wrong 
and/or if the shape of filter response function is incorrectly modeled since this will 
change the effective wavelength of the filters. We assume that this will induce fractional 
uncertainties in the effective wavelengths of [0.02,0.01,0.01,0.01,0.01] in [u, g, r, i, z] filters 
jsl. We note that this error will be a true systematic since the effect will be correlated 
between all source redshifts and colours. For comparison, for our QSO template, 
effective wavelengths will vary with ~ 2 — 3 % when varying the source redshift and 
dust parameters. Moreover, we assume a fractional error of 10 % in the dust map values 

§ A possible caveat not taken into account in this paper is that the dust map may be contaminated 
by the extragalactic background at low values (B. Menard, private communication). 
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Figure 2. Degeneracy between parameters p and Rv for SDSS colours u — g (dark 
red), g — r (red), r — i (orange) and i — z (yellow). The degeneracy lines correspond to 
a reddening equal to that of having p — 0.02 and Ry — 3.0. It is evident that multiple 
colours allow the degeneracy to be broken. 
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Figure 3. Dust extinction curves for different values of Ry shown together with the 
(arbitrary normalised) SDSS filter transmission curves. Left: Extinction curve for the 
CCM parametrisation. Right: Extinction curve for the Fitzpatrick parametrisation. 



We incorporate uncertainties using a Monte Carlo based approach as follows: 

• Pick a random value for the shift in effective wavelength for each filter and compute 
the corresponding colour offset. Also randomise the value for at each sky 
position within the given uncertainties. 

• Derive the likelihood over a grid of fitted parameters. In doing this, use the 
statistical sources of error from observational uncertainties and the intrinsic colour 
dispersion of the sources. 

• Repeat N times (we have used ~ 50) to obtain likelihood functions. 

• Compute the mean of the A^ likelihood functions. 

In this way, we take full account of the different correlations between sources of error, and 
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Figure 4. Left: Angular positions of the SDSS QSOs used in this analysis with 
darker colours corresponding to higher observed surface densities of QSOs, reflecting 
the given survey strategy. Right: Restframe QSO template spectrum shown together 
with arbitrary normalised SDSS filter transmission curves. 



the resulting confidence contours correctly incorporates both statistical and systematic 
effects. As noted in Sec. [H the dominating source of error is the possibility of shifts in the 
effective filter wavelengths, originating in uncertainties in the filter response functions 
and the source spectral templates. 



4. Quasar analysis 



It is well known that QSOs form a relatively homogeneous group of objects in terms 
of colours and spectral features. They also have large flux at short wavelengths which 
means that we can use blue colours, e.g., u — r also for high redshift sources, allowing 
us to effectively break the degeneracy between p and Ry as shown in Fig. [2j We use 
the fifth edition of the Sloan Digital Sky Survey (SDSS) Quasar Catalog [13f| consisting 
of 77429 objects. The area covered by the catalogue is 5 740 deg^, see Fig. HI 

QSOs have also been used successfully to constrain the amount and properties of 
dust in galaxies at cosmological distances p, B] and for intergalactic dust . 

Since we are primarily interested in the colours of QSOs, it is important to 
understand how the objects in the SDSS QSO catalogue are selected in order to assess 
whether there are selection effects that need to be taken into account. The QSO 
candidate selection is based on the position of sources in multidimensional SDSS colour 
space. In terms of colours, QSOs are identified as being outside the four-dimensional 



stellar loci points in [m, g^ r, i, z\ space plus a number of extra inclusion and exclusion 
regions. Since the selection is made after correcting for the extinction maps in Schlegel 
et al. 1^, assuming that these provide a fair model for the extinction, there is no reason 
that the selection should be severely affected by Milky Way dust extinction. 

Only QSOs in the redshift range 0.5 < z < 3.0 are considered, as QSOs outside this 
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Figure 5. Observed g — r colour as a function of galactic dust column density for 
SDSS QSOs at three different redshifts. Note the close to, but not perfectly, linear 
relationship between the observed colour and the dust density. 



range are less accurately described by the _&pectral template used, created by combining 



the HST radio-quiet composite spectrum [ij] with the SDSS median composite spectrum 



15[ as described in Ostman et al. [7|,|8[. Generally, all colours are relatively homogeneous 
at z < 2. Typically, at 2; ~ 1, we have a{u — g) = 0.15, a{g — r) = 0.13, a{r — i) = 0.093 
and a{i — z) = 0.11. However, when the Lyman-a emission is being redshifted 
into the observed filters, the dispersion increases. We thus cut the w-band data at 
z > 3000/1210 - 1 ~ 1.5 and the ^-band data at z > 4000/1210 - 1 ~ 2.3. 

In order to display the close to linear relationship between the observed QSO colour 
and the dust density, in Fig. [5] we show the observed QSO g — r colour as a function of 
galactic dust column density at redshifts z ~ 0.95, z ~ 1.72 and z ~ 2.48. 

Results for the full fit of p and Rv are shown in Fig. [6] for the COM dust 
extinction parametrisation (left panel) and the Fitzpatrick parametrisation (right 
panel). Solid colours correspond 95.4% CL. for two parameters. Our constraints are 
p = 0.0205 ± 0.0015 and Ry = 3.10 ± 0.43 (COM) and p = 0.0190 ± 0.0015 and 
Rv = 3.25 ± 0.37 (FTZ) (95.4% CL. for one parameter). Although consistent at the 
95.4 % C.L., compared to the r — i and i — z colours, we have on average less extinction 
in u — g and more extinction in g — r for the CCM parametrisation. For Fitzpatrick 
dust, compared to g — r and r — i, u — g has less extinction and i — z more extinction 
on average. 

Albeit with very limited resolution, multiple colour data allow us to also derive the 
dust extinction curve directly from the data, without the need to assume any specific 
dust extinction parametrisation. For each colour, we fit the slope, kxY, in Eq- El apply 
the necessary /^-corrections and normalise the extinction to be zero in the V"-band. In 
Fig. m we show the dust law as derived from QSO colours compared to the extinction 
curves for some fiducial dust models. Assuming the errors to be accurately represented 
by the uncertainty of the slope in the linear fit, the error bars are too small to be 
visible. However, this assumption is valid only in the case that the relation between the 
reddening and the dust column density indeed is linear and constant over the field. Also, 
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Figure 6. Constraints on p and Rv as derived from the observed colours of SDSS 
QSOs and the measured dust column density at the angular position of each QSO. 
Solid colours correspond 95.4% C.L. for two parameters. Left: Results for CCM 
dust corresponding to 95.4% C.L. for one parameter oi p — 0.0205 ± 0.0015 and 
Rv = 3.10 ± 0.43. Right: Results for Fitzpatrick dust corresponding to 95.4% C.L. 
for one parameter o{ p = 0.0190 ± 0.0015 and Ry = 3.25 ± 0.37. 
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Figure 7. Empirical dust law as derived from QSO colours (diamonds) compared 
to the extinction curves of CCM dust with p — 0.01, Ry = 5.0 (dotted curve), 
p = 0.02, i?y = 3.0 (dashed curve) and p = 0.03, i?y = 1.0 (dash-dotted curve). 



the dominating source of error of possible systematic shifts in the effective central filter 
wavelengths have not been included. Note that since the normalisation is arbitrary, only 
the shape of the derived extinction curve is relevant. 
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Figure 8. Left: Angular positions of the SDSS BCGs used in tliis analysis with 
darker colours corresponding to higher observed surface densities of BCGs. Right: 
Restframe BCG template spectrum shown together with arbitrary normalised SDSS 
filter transmission curves. Note that the same template is used for the LRG analysis. 



5. Brightest Central Galaxies analysis 

The brightest galaxies at the centre of galaxy clusters are empirically known be relatively 
homogeneous in colour. Since high luminosity early type galaxies are generally redder 
than low luminosity ones (usually attributed to higher metallicity), the colour dispersion 
can be further reduced using colour-magnitude relations. However, in this paper, we 
defer from employing any such calibration. 

We have used data from the MaxBCG catalogue of 13 823 galaxy clusters from 
the Sloan Digital Sky Survey Q. In order to extract the full photometric data for the 
brightest central galaxies (BCGs), we retrieved all objects within 6 arc sec from the BCG 
positions. We managed to identify and retrieve [u,g,r,i,z] data for 13 796 BCGs with 
redshifts 0.046 < z < 0.37. The angular positions of the BCGs are depicted in Fig. |H1 
Generally, all colours involving the u-band are not very well-determined. Typically (at 
z ~ 0.2), we have a{u-g) = 0.32, a{g-r) = 0.059, a{r-i) = 0.029 and a{i-z) = 0.042. 

We cut all colours involving the the u-band at z > 0.15 since the colour dispersion 
becomes very large at those redshifts. Results for fitting p and Ry to the BCG data are 
shown in Fig. M Our constraints are p = 0.0205 ± 0.0030 and Ry = 2.89 ± 0.88 (CCM) 
and p = 0.0177 ± 0.0020 and Ry = 3.42 ± 0.48 (FTZ) (95.4% C.L. for one parameter). 
Similar to the QSO constraints, compared to r — i and i — z that gives p = 0.02 for 
Ry = 3, we have less extinction in u — g and more extinction in g — r for the CCM 
parametrisation. For Fitzpatrick dust, g — r and r — i gives p = 0.02 at Ry = 3 while 
u — g has less extinction and i — z more extinction. In Fig. [lOl we show the empirically 
derived dust law. 
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Figure 9. Constraints on p and Rv as derived from the observed colours of SDSS 
BCGs and the measured dust cohimn density at the angular position of each BCG. 
Solid colours correspond 95.4% C.L. for two parameters. Left: Results for CCM 
dust corresponding to 95.4% C.L. for one parameter of p = 0.0205 ± 0.0030 and 
Rv = 2.89 ± 0.88. Right: Results for Fitzpatrick dust corresponding to 95.4% C.L. 
for one parameter of and p = 0.0177 ± 0.0020 and Ry = 3.42 ± 0.48. 
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Figure 10. Empirical dust law as derived from BCG colours (diamonds) compared 
to the extinction curves of CCM dust with p = 0.01, i?y = 5.0 (dotted curve), 
p = 0.02, Rv = 3.0 (dashed curve) and p — 0.03, Ry = 1.0 (dash-dotted curve). 
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Figure 11. Constraints on p and Ry as derived from the observed colours of 
SDSS LRGs and the measured dust column density at the angular position of each 
LRG. Solid colours correspond 95.4% C.L. for two parameters. Le/i; Results for CCM 
dust corresponding to 95.4% C.L. for one parameter of p = 0.026 ± 0.004 and 
Rv = 1-95 ± 0.35. Right: Results for Fitzpatrick dust corresponding to 95.4% C.L. 
for one parameter of and p = 0.019 ± 0.0025 and Ry = 3.16 ± 0.36. 



6. Large Red Galaxies 

Since Large Red Galaxies (LRGs) are also known to be homogeneous in terms of colours, 
we have extracted data for 114151 LRGs from the SDSS data base. Generally, all 
colours involving the u-band are not well-determined. Typically, we have (at z ~ 0.3): 
a{u - g) = 0.84, a{g - r) = 0.086, a{r - i) = 0.043 and a{i - z) = 0.056. We 
therefore do not use any colours involving the w-band in this analysis. Fitting for p 
and Rv gives the results shown in Fig. [TTl Our constraints are p = 0.026 ± 0.004 and 
Rv = 1.95 ± 0.35 (CCM) and p = 0.019 ± 0.0025 and Rv = 3.16 ± 0.36 (FTZ) (95.4 % 
C.L. for one parameter). Compared to the QSO and BCG analysis, the absence of 
the M-band pushes best fit values towards higher p and lower Rv- Assuming Rv = 3, 
compared to r — i and i — z that gives p = 0.02, we have more extinction in g — r for 
the CCM parametrisation. For Fitzpatrick dust, r — i gives p = 0.02 while g — r and 
i — z have more extinction. In Fig. [121 "we show the derived dust law. 

7. Discussion 

As seen in the previous sections, there are slight differences in the preferred values of 
p and Rv, depending on the dust parametrisation and the data set used. Generally, 
for the CCM parametrisation with Rv = 3, for r — i and i — z, we get p = 0.02. We 
have more extinction (higher p) in g — r and less in u — g. For Fitzpatrick dust, r — i 
gives p = 0.02 while g — r and i — z have more extinction and u — g less. How can we 
understand these differences? First, we need to understand the difference between the 
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Figure 12. Empirical dust law as derived from LRG colours (diamonds) compared 
to the extinction curves of CCM dust with p = 0.01, i?y — 5.0 (dotted curve), 
p = 0.02, Rv ~ 3.0 (dashed curve) and p — 0.03, Rv — 1-0 (dash-dotted curve). 



CCM and FTZ dust extinction parametrisations. In Fig. [131 we plot 

JWfTZ — fWcCM 
/(A) CCA/ 

The reason that the we need higher p for i — z for the FTZ model is that the extinction 
curve is flatter there compared to the CCM case. We also see that it is steeper for 
g — r which means that we get a lower p for that colour. That is, we need a dust 
law that is slightly steeper between for g — r and flatter for u — g compared to CCM. 
As for differences between the data sets, comparing BCGs and LRGs, we see that the 
differences can be explained by having a flatter template spectrum for LRGs in the 
region covering g — r, compared to the BCG template. Comparing BCGs and QSOs, 
differences can be resolved by either having the BCG template flatter for u — g and g — r, 
or the QSO template steeper in the same wavelength region. 

8. i^y- variations 

It is not unreasonable to assume that the distributions of grain sizes, and thus also the 
resulting value of Ry, is different in high and low density environments. Whether grains 
are generally larger or smaller in high density environments depend on whether grain- 
grain collisions are more likely to cause coalescence or shattering. We therefore check 
for any evidence for correlations between Ry and the dust column density by dividing 
our sample into subsamples with different column densities. Note however that a high 
column density not necessarily corresponds to a high physical density. In Fig. [HI we 
show constraints from dividing the QSO sample into two regions with low {D^ < 3 MJy 
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Figure 13. The difference between tlie CCM and FTZ dust extinction curves. 



sr ^) and high {D'^ > 3 MJy sr ^) dust column density. Our constraints for CCM are 
p = 0.020 ± 0.0016 and Ry = 3.06 ± 0.56 {D^ < 3 MJy sr'^) and p = 0.019 ± 0.0019 
and Rv = 3.09 ± 0.65 {D^ > 3 MJy sr-^). For Fitzpatrick dust, p = 0.019 ± 0.0016 and 
Rv = 3.28 ± 0.37 {D^ < 3 MJy sr^^) and p = 0.018 ± 0.0018 and Ry = 3.23 ± 0.41 
{D^ > 3 MJy sr~^) (95.4% C.L. for one parameter). This simple test does not provide 
evidence for any correlation between the dust column density and dust properties. 



8.1. cr{Rv) and a{A 



V) 



If dust has spatially varying properties, depending on the size of the regions, we would 
expect the colour dispersion to be bigger behind regions with high column density D^, 
in addition to being redder. For filters X and Y, we can write the colour E{X — F ) as a 
combination of the intrinsic colour Ei{X — Y) and the dust induced colour Efi{X — Y), 



E{X -Y) = E,{X -Y) + E,{X - Y) , 
and the corresponding colour dispersion 

a[E{X - Y)]' = a[E,{X - Y)]' + a[E,{X - Y)]' 



where 



and 



Ax 
Av 



Ea{X -Y)=[-^-^]Ay = fxY{Rv)A 



a[E,{X - Y)\' 



Ay 

Av 



Rv dfxY 



R 



V 



+ 



Av 



JxY dRv 

= g[X, Y, Rv, a{Rv),a{Av)])E,{X - Yf 
In Fig. [151 "we plot the variance a[E{X — Y)]"^ as a function of Ed{X 



(9) 



(10) 



(11) 



E,{X - Y)\12) 



Yy for each 



standard colour of the QSO sample, expecting the variance at zero Ed{X — Y)^ to give 
the intrinsic colour variance a[Ei{X — Y)]'^ and the slope to be related to the scatter 
in Rv and Av through the function g[X,Y, Rv,a{Rv),(j{Av)] (using the CCM dust 
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Figure 14. Constraints onp and Ry when dividing data into regions with low {D^ < 3 
MJy sr^^) and liigh {D^ > 3 MJy sr~^) dust column density. There is no evidence 
for any correlation between the dust column density and dust properties. Upper left: 
Results for CCM dust for all regions with < 3 MJy sr^^ corresponding to 95.4% 
C.L. for one parameter of p = 0.020±0.0016 and Ry — 3.06±0.56. Upper right: Results 
for CCM dust at D"^ > 3 MJy sr^^ giving p = 0.019 ± 0.0019 and Ry = 3.09 ± 0.65. 
Lower left: Resuhs for Fitzpatrick dust at I?"^ < 3 MJy sr'^ giving p = 0.019 ±0.0016 
and Ry ~ 3.28 ± 0.37. Lower right: Results for Fitzpatrick dust at > 3 MJy sr^^ 
giving p = 0.018 ± 0.0018 and Ry = 3.23 ± 0.41. 



extinction parametrisation). As can be seen, there are no obvious correlations between 
a[E{X — Y)]'^ and Ed{X — Yy (or equivalently the dust column density D^). Assuming 
ajiy/Rv and (Ta^/^v to be constant (although the results depicted in Fig. [15] indicate 
that this assumption may be too simplistic), we can constrain their values from the 
fitted straight line slopes. Results from such an analysis are shown in Fig. [161 The 
bluer colours show no increase in the colour variance with dust column density, whereas 
i — z show a slight increase in the variance indicating a possible variance in aAY/^v at 
the 30 % level. This said, there are a few caveats that should be remembered: First, this 
is under the assumption that the relation between the variance and the colour indeed 
is linear. Second, results for i — z are not consistent with the bluer colours that show 
no signs of a scatter in the dust properties (although it is possible that the scatter only 
manifests itself at longer wavelengths). Third, there is a (weak) binning dependence 
when calculating the dispersion as a function of Ed{X — Y)"^. To summarize, there is no 
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Figure 15. The variance a[E{X — Y)]'^ as a function of Ed{X — for each standard 
colour of the QSO sample. The total variance at zero Ed{X — Y)'^ is the intrinsic 
colour variance a[Ei{X — Y)]^ and the slope is related to the scatter in Ry and Ay- 



evidence in the QSO data for any additional scatter in the observed colours in regions 
with high column density, indicating that dust properties are relatively homogeneous 
averaged over the size of spatial regions probed in this paper, except perhaps for a hint 
of such an increased scatter at the reddest optical wavelengths probed by the QSO data. 

9. Summary 

By investigating the correlations between dust column density as inferred from infrared 
data and the observed colours of celestial objects at cosmological distances with small 
colour dispersion, we constrain the properties of Milky Way dust. Results from QSOs, 
BCGs and LRGs colours are broadly consistent, indicating a proportionality constant 
between the reddening E{B — V) = As — Ay and the dust column density 
(given in units of MJy/sr) of p = E{B — V)/D'^ ~ 0.02 and a reddening parameter 
Ry = Ay/E{B — V) ~ 3 with relative uncertainties of the order 10%. In spite of the 



increased statistics, the fractional error of p is similar to that of Schlegel et al. [12[, the 
reason being that we simultaneously fit for Ry and that the uncertainty is dominated 
by the possibility of systematic shifts in the effective filter wavelengths. 

The data does not provide any evidence for any dependence of the dust properties as 
a function of dust column density, neither for any sizable scatter in the dust properties on 
the scales probed by the data, except a hint of scatter a Ay /Ay > 0.3 in the extinction 
properties of dust at the longest wavelengths probed by the data (corresponding to 
i — z). We have also checked that within the uncertainties, there is no evidence of 
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Figure 16. Constraints on constant /Ry and (jav l^v derived from the linear fits 
depicted in Fig. [15] The coloured regions correspond to 95.4% CL (red) and 68.3% 
CL (yellow), both for two parameters. 



any systematic difference between the dust extinction properties in the northern and 
southern hemisphere. Assuming the dust maps of Schlegel et ah [l2|| to be correct, it is 
thus fair to assume that we can correct for Milky Way dust by applying a proportionality 
constant p = E{B — V)/D'^ ~ 0.02 with reddeni ng p arameter Ry ~ 3 to the dust maps, 
in accordance with the results of Schlegel et al. 12| . 
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